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ABSTRACT 

■ Pulsar emission should primarily come from the magnetic separatrix... Combining 

. theory and observations, we show that force-free electrodynamics (FFE) gives an 

, accurate description of the large-scale electromagnetic field in the magnetospheres 

I of Crab-like pulsars. A robust prediction of FFE is the existence and stability of a 

I singular current layer on the magnetic separatrix. We argue that most of the observed 

pulsar emission comes from this singular current layer. 
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■ 1. Introduction 
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' Two different things are done in this paper. First, we show that simple and pure theory - 

■ force-free electrodynamics (FFE) - provides an accurate quantitative description of the large-scale 
electromagnetic field in pulsar magnetospheres. Second, we use certain FFE results to propose 
that pulsar emission is mostly generated in the vicinity of the magnetic separatrix, which is the 

Q . surface separating the regions of closed and open magnetic field lines. 

H ; 

^ I FFE describes electromagnetic fields in plasmas, but the only degrees of freedom in FFE are 

the electromagnetic fields of a special force-free geometry. The plasma is described implicitly. 
FFE makes no predictions regarding plasma emission. 
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Stationary FFE equations describing pulsar magnetospheres have recently been solved 
(Contopoulos, Kazanas, Fendt 1999, Gruzinov 2005, 2006, Spitkovsky 2006). The shape of the 
pulsar magnetosphere, the singularity structure, and the spin-down power have been calculated (in 
a sense, see §3). But these results hang in the air; being pure FFE results, they make no directly 
testable predictions ^. 

We think that FFE has actually found the primary site of pulsar emission. We propose 
that pulsar emission is due to a small dissipation effect on top of the ideal dissipationless 
FFE. The dissipation is basically Joule heating. The current layer on the separatrix consists 
of counter-propagating beams of electrons and positrons of formally infinite densities. For any 
"reasonable" anomalous resistivity, the Joule heating in the singular current layer exceeds the 



^One may even say that FFE does make one prediction - the spin-down rate should scale as the cube of the spin. 
This is known to be wrong. But no satisfactory explanation of the anomalous breaking of pulsars exists today. We 
will not count this as an argument against FFE. 
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Joule heating in the bulk. Then the singular current layer should also dominate the pulsar 

emission. Our theory is obviously incomplete - we do not calculate the emitted power spectrum 
and the light curves. But the case for the magnetic separatrix being the primary source of pulsar 
emission looks strong. 



2. FFE 

Here we formulate FFE and discuss its meaning and conditions of applicability. Then we show 
that FFE provides a valid description of the large-scale electromagnetic field in the magnetospheres 
of Crab-like pulsars. 

In FFE, as in all plasma models, Maxwell equations are used for the electromagnetic field: 

5tB = -VxE, atE = VxB-j. (1) 

To close this system, one needs to calculate the plasma current density j. In standard plasma 
physics this is done by solving the equations of motion of charged particles. In FFE one postulates 
the following Ohm's law: 

. (B-VxB-E-Vx E)B + (V • E)E x B 

J = ^2 ■ (2) 

Equations (1,2) form a system of evolutionary nonlinear partial differential equations - FFE. FFE 
is Lorentz covariant. 

The FFE Ohm's law was designed (Gruzinov 1999) to maintain the force-freedom: 

pE+jxB = 0, p = V-E. (3) 

One can show that if the force-free constraint (3) is satisfied at the initial time (via E ■ B = at 
the initial time), the FFE time evolution will maintain the force-free constraint (with E • B = 0). 

Consider the realm of applicability of FFE. If the stress-energy tensor of the system (plasma 
plus electromagnetic field) is dominated by the electromagnetic part, the force-free constraint is 
automatically satisfied. But it is unclear in which conditions would the plasma arrange the fields 
into a forcc-frcc configuration instead of accelerating the charges up to the energy equipartition 
with the field. 

The force-free constraint takes a simple form in particular inertial frames. For arbitrary 
electric and magnetic fields, E and B, at a given spacetime event, one can always find inertial 
frames with E parallel to B at this event. In such frames, the force-free condition means that 
(i) the current j is along B, (ii) the electric field E vanishes. We will use this formulation of the 
force-freedom in our discussion of pulsar magnetospheres. 

Consider a pulsar with the following Crab-like nominal parameters. Radius of the star 
Ro = 10 km, angular velocity J7 = 200 s~^, spin-down power L = 5 x 10^^ erg/s. The nominal 
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magnetic field at the pole Bq = 5.8 x 10^^ G (for axisymmetric pulsar, L fi^fl^/c^, fi is the 
dipole). The light cylinder radius is Ri^ = c/Q. = 1500 km. The nominal magnetic field at the 
light cylinder B\c = i?o(-Ric/-Ro) = 1-7 ^ 10^ G. The nominal particle density near the star is 
no = rii?o/(27rce) = 1.3 x 10^^ cm^^ (Goldreich, Julian 1969, derived from the typical charge 
density required by stationary FFE, §3). The nominal particle density near the light cylinder is 
nic = 3.8 X 10^ cm-3. 

The first force-free condition - current flows along B - must be satisfied. Even near the light 
cylinder, the cyclotron cooling time of electrons and positrons, (3mc)/(4rgS[^) = 9 x 10~^ s {r^ 
is the classical electron radius), is much smaller than the pulsar period, 27r/J7 = 3.1 x 10~^ s. 
Most of the charges must be in the lowest Landau state, with the residual perpendicular current 
corresponding to magnetization ~ nicf^B = 3.5 x 10~^^ G, fiB is the Bohr magneton. We see that 
for any possible multiplicity (the ratio of true to nominal particle density), the perpendicular 
current is negligibly small. 

Actually the leading cause of current deviation from the magnetic field lines is the curvature 
of the magnetic field lines ^. To estimate this effect, one must calculate the ratio of the centrifugal 
pressure of the outflowing plasma to the magnetic pressure. Near the light cylinder, this reduces 
to the ratio of energy densities. The nominal value of this ratio, calculated assuming mildly 
relativistic outflow with small multiplicity, is ~ nicmc /Bl = 1.1 X 10-^2 This must remain small 
for any reasonable multiplicity and Lorentz factor of the outflowing plasma, because the plasma 
outflow with high Lorentz factor would produce more curvature radiation than what is actually 
observed in the Crab pulsar (see below). 

The second force-free condition - zero electric field along the magnetic field - is more 
stringent. One can give a simple theoretical proof that the parallel electric field vanishes for radii 
^ 1000 km. But this is not good enough to claim that FFE is quantitatively correct. We need 
FFE to work all the way to few light cylinder radii. Even though we do not know how exactly 
it happens, we will show, using the observational upper bounds on TeV Crab emission, that the 
parallel electric field must nearly vanish near the light cylinder too. 

First consider the ^ 1000 km zone. Here a clear-cut mechanism is available that would screen 
out the parallel electric field (Sturrock 1971, Rudcrman, Sutherland 1975): (i) parallel electric 
field accelerates electrons and positrons, (ii) accelerated charges emit curvature photons, (iii) 
curvature photons produce electron-positron pairs on the background magnetic field, (i) parallel 
electric field accelerates electrons and positrons... - leading to an avalanche that screens out the 
parallel electric field. We do not claim that this mechanism does operate in real pulsars. We only 
point out that this mechanism kills the parallel electric field, if the field is not actually killed by 
some other more powerful mechanism. 

For the numerical estimate, suppose that FFE does not apply starting from radii ^ R. Then 



thank Jonathan Arons and Anatoly Spitkovsky for explaining this to me. 
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the parallel electric field at radius R is close to the characteristic FFE value, E ^ {QR/c)B, where 
the characteristic magnetic field at radius i? is -B ~ Bo{R/Rq)~^. This electric field will accelerate 
electrons and positrons to such Lorentz factors 7, that the curvature radiation power balances 
the electric power, e^cy ~ ceE, giving the terminal Lorentz factor 7 ~ (27rnoi?g)^/^ = 10^. 
Electrons and positrons moving at such Lorentz factors emit curvature photons of energy 
e ~ 'y^ch/R. These photons will produce electron-positron pairs on the magnetic field B with the 
mean free path ^ R, provided x = {eB)/{2mc^Bq) > 0.07, where Bg = vn?c^/{eh) = 4.4 x lO^^G. 
This gives the maximal radius for the pair production avalanche R < 780 km. 

Now consider the zone around the light cylinder. Again assume that FFE does not apply. 
Then the parallel electric field will accelerate all available charges to Lorentz factors 7 = 10^. 
These charges will emit curvature photons of energy eic ~ ^'^ch/R\f. = 0.13 TcV. Since we know 
that FFE does apply just inside the light cylinder, the density of the charges near the light cylinder 
cannot be much smaller than ^ nic, and therefore the pulsed Crab emission at ^ 0.1 TeV cannot 
be much smaller than ~ R^^n\ceB\cC L = 5 x 10^^ erg/s. This is in gross contradiction with 
observations. The upper bound on Crab's pulsed emission at 100 Gev is about ~ 3 x 10^^ erg/s, 
assuming 2 kpc distance (Oser et al 2001, extrapolation of Aharonian et al 2004). 

In summary: the current flows along the magnetic field, the parallel electric field vanishes. 
This means that the plasma is force-free, and the dynamics of the electromagnetic field reduces to 
FFE. 



3. Pulsar FFE and Pulsar Emission 

In covariant form, the force-free constraint (with E • B = 0) gives F^^^ jy = 0, where F 
is the electromagnetic field tensor, and j is the 4-current. It follows that the electromagnetic 
stress-energy tensor is conserved, dyT^^^ = 0, and FFE conserves energy £, momentum P, and 
angular momentum M 

£ = j d^r ^ , P = J (fr ExB, M = j d^r r x (E x B). (4) 

The stationary (that is co-rotating) FFE magnetosphere of the pulsar of angular velocity fl is the 
stationary point of free energy (Gruzinov 2006) 

£' = £ -n-M = J d^r{ ^"^^ - n • r X (E X B) ) (5) 

under arbitrary variations JE of the electric field and iso-topological variations 5B of the magnetic 
field 

5B = V X ((5C X B). (6) 

Performing the variations, one gets the pulsar magnetosphere equation: 

B X V X (B + V X (V X B)) = 0. (7) 
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Here V = f2 x r is the co-rotation velocity, and the electric field is E = —V x B ^. The 

magnetosphere rotates at angular velocity fl. The actual time evolution probably does take the 
magnetosphere into this "minimum energy" state (Spitkovsky 2006). 

The fact that the stationary magnetosphere equation (7) follows form the iso-topological 
variational principle has two important consequences. 

Stationary magnetospheres, that is solutions of (7), exist for any prescribed topology of the 
magnetic field lines. Pure FFE cannot choose among these topologies. In particular, one can have 
arbitrary current flowing in the closcd-ficId zone. True, this current must dissipate, but it might 
be driven by the kinematic processes which produce the plasma. FFE cannot predict the pulsar 
power and magnetosphere, until the plasma production mechanism is understood quantitatively. 

Another consequence of the iso-topological minimization is the existence and stability of 
the singular current /charge layer along the magnetic separatrix (Gruzinov 2006). The singular 
current is known to exist in the axisymmetric magnetosphere (Contopoulos, Kazanas, Fendt 1999, 
Gruzinov 2005). The singular current layer is also seen in the generic inclined magnetosphere 
(Spitkovsky 2006). 

For the axisymmetric pulsar, assuming pure dipole on the surface of the star and assuming 
zero current in the closed zone, one can calculate the surface current and charge on the magnetic 
separatrix numerically with about 10% accuracy. For simplicity, consider only the region well 
within the light cylinder, r <^ 1, z <^ 1, where r is the cylindrical radius and z is the vertical 
cylindrical coordinate; we use natural "pulsar units" {Q = = c = I). 

The magnetic separatrix is a surface given by r'^/z^ ~ 1.25. The Goldreich-Julian charge 
density inside the separatrix is p = 2B = 4/z^. The current density on the magnetic axis (the 
line r = 0) is J !=a 3.8/ z^, comparable to to p. The current density decreases for growing r; the 
current reverses near the separatrix, but the total reverse current flowing inside the separatrix 
is small. There is no poloidal current outside the separatrix. Most of the reverse current flows 
in the singular current layer on the separatrix. The surface density of the separatrix current is 
i^l/r^ (0.8r)/z3. 

The surface charge of the separatrix, a ~ 0.2r, is much smaller than the surface current. This 
means that the separatrix current is formed by the counter-propagating surface flows of electrons 
and positrons. 

We will assume that the pulsar emission is primarily due to the anomalous Joule heating, and 
compare the Joule power of the separatrix and the Joule power of the bulk inside the separatrix. 
Since FFE was shown to be a good approximation, the width of the current layer on the separatrix 
6 due to kinetic non-FFE processes should remain small, (5 <C r. The current density in the 
separatrix layer jg ~ {'r'/S)ji, » ji,, where ji, is the current density in the bulk. The ratio of the 



^Which gives the Goldreich-Julian charge density used in §2. 
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Joule heating in the bulk to the Joule heating in the current layer 

Qb Ebjbv'^ El, 
Qs EsjsvS Eg ' 



(8) 



where E denote parallel electric fields in the bulk and in the separatrix current layer. For any 
"reasonable" anomalous resistivity, the electric field in the separatrix current layer should be much 
greater than the electric field in the bulk, because jg » jh- Then 

— ~ — < 1. (9) 

Given that the Joule heating in the bulk is much smaller than the Joule heating in the separatrix 
current layer, we propose that the pulsar emission is dominated by the separatrix current layer. 
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